Astronomy & Astrophysics manuscript no. paper 


© ESO 2008 


February 5, 2008 





A Method for Detection of Structure* 

M.Gustafsson 1 , J.L.Lemaire 2 , and D.Field 1 



O 
O 

(N 

C 



> 
o 
in 

o 

O 



1 Department of Physics and Astronomy, University of Aarhus, DK-8000 Aarhus C, Denmark 
e-mail: maikeng@phys . au . dk 

2 Observatoire de Paris & Universite de Cergy-Pontoise, LERMA & UMR 81 12 du CNRS, 92195 Meudon, France 

For Main Journal, Diffuse Matter in Space. Received: Accepted: 

ABSTRACT 

Context. In order to understand the evolution of molecular clouds it is important to identify the departures from self-similarity associated with 
the scales of self-gravity and the driving of turbulence. 

Aims. A method is described based on structure functions for determining whether a region of gas, such as a molecular cloud, is fractal or 
contains structure with characteristic scale sizes. 

Methods. Using artificial data containing structure it is shown that derivatives of higher order structure functions provide a powerful way 
to detect the presence of characteristic scales should any be present and to estimate the size of such structures. The method is applied to 
observations of hot H 2 in the Kleinman-Low nebula, north of the Trapezium stars in the Orion Molecular Cloud, including both brightness and 
velocity data. The method is compared with other techniques such as Fourier transform and histogram techniques. 

Results. It is found that the density structure, represented by Hi emission brightness in the K-band (2-2.5/mi), exhibits mean characteristic 
sizes of 110, 550, 1700 and 2700AU. The velocity data show the presence of structure at 140, 1500 and 3500AU. Compared with other 
techniques such as Fourier transform or histogram, the method appears both more sensitive to characteristic scales and easier to interpret. 
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1. Introduction 

The velocity and density fields of molecular clouds have been 
observed to show very complex structure. This property sug- 
gests that structure is dominated by turbulence and this has 
been verified by both observatio n s and numerical simula tions 
(e.g. iFalgarone &Phliii^ fl99(it iRoldvrev et alJ 
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also the extensive re views by Elmegreen & Scalo (2004) and 
Mac Low & Klessen (2004). Clouds dominated by turbulence 
lack any characteristic scale and are inherently self-similar 
with a variety of fractal dimensions. This has been observed 
over a broad range of scales, 0.02pc-100pc (e.g.lLarsonl ll 981; 
IChappell & Scalol200lUOssenkopf & Mac Lowl2002l) . 

Self-similarity and lack of characteristic scale in molec- 
ular clouds must however break down at some dimension or 
range of di mensions a ssociated with star formation, as demon- 
strated in IVannier et all {2001) and Ipustafsso n et alJ J2OO6). 
Turbulence creates local clumps of gas and this process, known 
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as turbulent fragmentation, is considered the root cause of ob- 
served density structure and may give rise to large density con- 
trasts. Star formation begins to take place at smaller scales in 
gravitationally unstable clumps. The distribution of sizes of 
these clumps is directly related to the initial mas s function 
(IMF) for star formation (Padoan & Nord lumfeool . The pre- 
dominantly turbulent era, characterized by lack of preferred 
scale, is expected to end when gravitation becomes dominant 
in clumps. At this stage, the density and velocity fields are 
not fractal, but will contain characteristic scales. Hence certain 
scales may be over-populated or under-populated relati ve to 
the model of a simp le Kolmogorov-type cascade ( iKolmogorovl 
ll94lUFrischlll995l) . The range of scales, at which clumps be- 
come gravitationally bound and turbulence looses its kinemat- 
ically dominating effect, is important for both theories of star 
formation and of turbulence. 

Characteristic scales may also be imposed on the medium 
by other mechanisms than gravitational collapse. In particular 
energy may be injected into the medium at the largest scales by 
supernova explosions or at the tens of thousands of AU scale 
by outflows from massive stars or from low mass stars at scales 
of a few thousand AU. In order to understand the evolution of 
molecular clouds it is essential to identify the presence of these 
energy injection processes and the scales at which they occur. 
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The purpose of the present paper is to introduce a tech- 
nique, based on structure functions (SFs), to identify charac- 
teristic scales in the interstellar medium - or indeed in any 
medium. The technique appears more sensitive to characteristic 
scales, for example in spatial brightness or velocity, than meth- 
ods presently used. The method is relatively easy to interpret 
in terms of preferred scale sizes and is as easy to implement as 
existing methods. 

Experience has proven that it is non-trivial to establish 
the scale(s) at which self-similarity begins to break down and 
an array of techniques as broad as possible should be made 
available. For exam ple, four different methods were used in 
IVannier et alJJ20Qll) in order to demonstrate that self-similarity 
breaks down at star-forming scales of around 500 - 1000 AU 
in a small portion of the Orion Molecular Cloud (OMC1). 
The present contribution should be seen in this light, not as 
a panacea for structure identification but rather as another tool 
for this purpose. 

The methods curren tly available are those of 
(i) iBlitz & Williams! dl997h ." who introduced a method of 
degrading the spatial resolution and establishing the pres- 
ence of preferred scales from histogram s of the de graded 
images, a met hod u s ed for example in IVannier et all ll200ll) 
and Lacomb e et alJ (|2004), (ii) Fourier transform power 
spectra, used to pick out f r equen cies and hence scales in 
images (e.g. | Va nnier et all 1200 lb . (iii) clump decompo- 
sition Jstutzki & Gueste nl Il990h . (iv) A-variance analysis 
(Stutzki et al. 1998 ; iBensch et alJEoOll) . (v) wavelet transfor- 
mation (e.g. Farge 1 992) . 

The non-fractal nature of molecular clouds at the star- 
forming scal e of ~ 500-1000 AU was first detected in 
IVannier et alJ J200ll) for Orion. Essentially the same pre- 
ferred scales b e tween 500 and 1100 AU are also reported in 
lLacombe et alJ J2004I) . Non-fractal nature has al so been re- 
ported at larger scales in Blitz & Williami il 19971) who found 
a preferred scale of 0.25-0.5 pc (5-10xl0 4 AU) in Taurus. 
Material may however be self-similar down to a few x 10~ 4 pc 
(<100 AU) in diffuse material, with clumps of mass ~ 1O~ 7 M , 
as found in the photodisso ciation region around the B-star 
HD37903 in NGC2023 kouan et alJl997l) . Here we show that 
structure functions (SFs) or local derivatives of SFs provide 
a powerful way to detect preferred scales in a medium. The 
method is tested on observations of H2 emission in OMC1. 

2. Method 

The SFs of order p of a spatially resolved parameter A are de- 
fined as 

S p (r) = <|A(r') - A(r' - g)\ p ) = <\AAf) (1) 

where the average is taken over all map positions r' and all dis- 
tances 5 such that \8\ - r. The parameter A can for example be 
velocity or brightness. The SFs are measures of the spatial cor- 
relations of the parameter A. A completely random distribution 
would result in a constant structure function as a function of r, 
whereas for example a gradient in the image would lead to a 
structure function increasing with r. 



Structure functions have traditionally been used in turbu- 
lence theory to quantify the properties of the turbulence. In a 
predominantly turbulent medium where the velocity and den- 
sity fields are self-similar or fractal, the SFs are well described 
by power laws, S p (r) oc r*^ in the inertial range between 
the d riving and dissipat ion scale. This is foun d from both the- 
ory (Kolmogorov 1941; Boldvrev et al. 2002b) and simulations 
(e.g. lBoldvrev et alJ2002aHKritsuk & NormarJ2004 and sim- 
ply conveys that in turbulent regions, the velocity field, say, is 
spatially correlated in such a way that larger separations corre- 
spond to larger velocity differences on average. This is in fact 
a manifestation of the turbulent energy cascade originally pro- 
posed by Kolmogorov. 

On the other hand, if the region in question is not self- 
similar, but contains structure with a preferred scale, the SFs 
will deviate from power-laws. It is shown here that these devi- 
ations are a sensitive test of the presence of structure and may 
be used to reveal the size of structures. If the region contains 
only a single structure in the map of parameter A of a certain 
size, R, it is evident that correlations will only persist up to 
scale R. In the structure function this will show up as an in- 
crease of S p (r) at scales < R and a constant value at scales > 
R. The same behaviour will be found if there is a large num- 
ber of structures with the same size randomly distributed in the 
map. When clumps with two or more distinct preferred sizes 
are present, which is most likely the case in star forming re- 
gions, the structure function will show this effect for all sizes 
involved. The effect on the SFs will become successively less 
pronounced at any specific size as more scales are identified. In 
a real system, structures will most likely tend to group around a 
preferred scale with some characteristic deviation from a mean 
value. This aspect turns out to be more of a gloss on the present 
description than a fundamental point but can lead to system- 
atic errors in structure size determination. This is described in 
Sect.l4~4l 

To test the ideas set out above, three maps have been con- 
structed with i) a single circular clump with a diameter of 40 
pixels and an intensity distribution given by the paraboloid 
function 

tt \ rn (x ~ x ) 2 + (y - y Q ) 2 

f(x,y) = C(l j ) (2) 

where r is the radius of the clump, (xo,yo) is the centre position 
and C is the peak intensity, ii) 100 randomly distributed clumps 
with the same characteristics as in i) and iii) 100 clumps with 
a diameter of 20 pixels and 20 clumps with a diameter of 100 
pixels where again all clumps were randomly distributed and of 
paraboloid shape. In all cases the clumps have been distributed 
on a square 1024 by 1024 grid. Similar maps were also con- 
structed using conic shapes (f(x,y) = C(l — — • t °^ +( -L 2£l)) 
for the intensity distribution in the clumps. 

The SFs of order 1-5 for case (i) are shown in Fig.^. The 
SFs become roughly constant at r > 30 pixels close to the ex- 
pected value of 40. The SFs for case (ii) are essentially the 
same. They are not shown here. The effects on the SFs caused 
by the presence of clumps are more clearly seen at higher or- 
ders, p. 
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Fig. 1. a) Structure functions of order 1-5 (labelled accordingly) for 
a simulated map, case (i) (see text), b) Local logarithmic derivatives 
of the SFs shown in a) Inset: Blow-up of r=29-45 pixels, c) Local 
logarithmic derivatives of SFs of order 1-5 derived from the simulated 
map of case (ii). 



Preferred scale sizes in a map or an image result in changes 
in the logarithmic slope of SFs. Therefore the local logarithmic 
derivatives, dlogS p (r)/dlogr, highlight the deviations from 
fractality and thus we propose to use these derivatives to de- 
tect such scales. As seen above the presence of preferred scales 
causes a decrease in the slope of the SFs making it roughly 
constant near the relevant size. This results in a plateau in the 
logarithmic derivatives if larger scales are not present, or a min- 
imum if such scales are present. The former case is demon- 
strated in Fig. [fl? for case (i) above, that is, for a single clump, 
and in Fig.QJ for case (ii), that is, for 100 such clumps. A min- 
imum is seen at ~ 35 pixels in both Figs, \I]p and c, close to 
the imposed structure size of 40 pixels (see above). The inset 
to Fig. \I]p displays the derivatives around r~ 35 in detail and 
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Fig. 2. Local logarithmic derivatives of SFs of order 1-5 from 
the simulated map of case (iii) (see text). Inset: the logarithm 
of the 5th order structure function of case (iii). 

shows that the position of the minimum moves towards a lower 
value of r when p is increased. This occurs because higher or- 
ders of the structure function give more weight to higher val- 
ues of AA (see Eq. ([0), thus depressing the outer regions of 
clumps, which is the location of low values of AA in these sim- 
ulations. Thus estimates of scale sizes from high order SFs rep- 
resent lower bounds to the scale sizes and the estimated values 
can be 10-20% too low. In Fig.^;, in addition to the first plateau 
or weak minimum at ~ 35 pixels, there is a secondary plateau 
at ~ 75 pixels, indicating the presence of a scale size approx- 
imately twice as large as the clumps in the simulation. This is 
most likely due to overlapping of the randomly placed indi- 
vidual clumps. By implication this method cannot distinguish 
between one large clump or several smaller clumps nearly co- 
inciding in the same line-of-sight. 

The local logarithmic derivatives for p=l-5 are shown in 
Fig- HI for case (iii) involving 100 clumps of diameter 20 pix- 
els and 20 clumps of diameter 100 pixels. Here a local mini- 
mum is seen at ~ 17 pixels and ~ 100 pixels, representing the 
scales of 20 and 100 pixels included in this test image. Note 
that two scale sizes will only show clearly if the scales are well 
separated. If they are too close the larger scale will dominate 
and suppress the smaller. In real data this means that a de- 
tected scale size might conceal other smaller scales (see also 
Sect. I4.4t . The inset to Fig.[2]shows the corresponding struc- 
ture function for 5th order, showing plateaux at ~ 15 pixels and 
at >~80 pixels. 

The same tests were performed on the maps of clumps with 
conical shaped intensity distributions and yielded similar re- 
sults to those above. This suggests that effects seen in the struc- 
ture functions are independent of the shape of the clumps. This 
is of course only valid for clump geometries which can at least 
crudely be characterized by a single dimension. Tests were also 
carried out on clumps with an ellipsoid shape with randomly 
oriented axes. These are characterized by two sizes, that is, 
those of the major and minor axes and the results are very sim- 
ilar to those of clumps with two different sets of diameters, that 
is case (iii) above. 

An additional property is that in the presence of real struc- 
ture the derivatives of structure functions of all orders tend 
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to congregate around the same value at a local minimum or 
plateau. This provides a further characteristic which aids the 
detection of structu re. In numerical si mulations of fully devel- 
oped turbulence, Biferale et al. ( 2004) noted a similar tendency 
for the local derivatives to accumulate at a certain value. 

A further point in assigning the presence of structure is the 
question of how large a structure may be detected for a given 
size of image, so called "edge-effects". In case (iii), in Fig. 2, 
there is already evidence for apparent structure at ~ 250 pixels 
which is not in the data. Edge-effects have been studied here us- 
ing fractional Brownian motion structures (fBm-fractals). We 
have calculated SFs and logarithmic derivatives for a number 
of such fractal images. For these, the SFs are expected to be 
pure power laws and the derivatives should be independent of 
r. The value of the logarithmic derivative depends on the choice 
of power spectral index and the order of the structure function 
involved, being proportional to the latter. In passing we note 
that fractal systems cannot therefore yield logarithmic deriva- 
tives of SFs which congregate at a particular value. 

An example of these test calculations is seen in Fig.[3]where 
SFs of order 1-5 and logarithmic derivatives for a fractal image 
of size 2048x2048 pixels are shown. The SFs may be approx- 
imated by power laws save at the largest scales. The logarith- 
mic derivatives are nearly constant up to r ~ 700 pixels, any 
deviations from constancy below this value being due to the 
imperfect fractal nature of the simulations arising from pixe- 
lation. At around 700 pixels, a local minimum is found and 
the derivatives for p=3-5 congregate around the value of 0.4, 
apparently indicating a preferred structure size. The same anal- 
ysis has been carried out on fractal images of size 1024x1024, 
512x512 and 256x256 pixels. In all images apparent structure 
at scales of ~ 1/5 - 1/3 of the size of the map was found. Hence 
increasing the size of the map also increased the absolute value 
of the apparent structure size. Thus we conclude that the ap- 
parent structure in the fractal images is an edge-effect and that 
detection of scale sizes much greater than 1/5 of the size of the 
map should not be attempted as such scales cannot be reliably 
identified. 

3. Real data for test purposes 

In order to illustrate the method of structure analysis described 
here with real data, we use results for the emission from vibra- 
tionally excited H2 in OMC1, where excitation occurs largely 
through shock excitatio n (with a small photon excitation com- 
ponent: see e.g. Kristensen et al. 2003 1). The relevant d ata ar e 
reported in iG ustafsso n et alJ d200 3h. ILacombe et alJ (|2004), 
iNissen et all d2005l) and Kristen sen et al.1 d2006l) . These three 
datasets refer only to the hot component of the gas, above 
~ 1000K, but nevertheless form valuable datasets for the iden- 
tification of structu re since they probe v ery small scales, from 
as low as 55 AU (ILacom be et al. 2004), and also involve the 
major cons tituent of the gas, albeit in an excited state. 

Data in Gusta fsson et al.ld2003l) were obtained at the CFHT 
with a Fabry-Perot interferometer in conjunction with adaptive 
optics (GriF, Clen et et al.l 120 02) and contain both brightness 
and velocity information. Observations were performed in the 
NIR K-band by scanning the v=l-0 S(l) H2 emission line at 
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Fig. 3. a) Structure functions of order 1-5 (labelled) of fractal image, 
b) Corresponding local logarithmic derivatives of structure functions 
of order 1-5. 



2.121/im. The field of view was 36"x 36"and the pixel scale is 
07035 (1024x1 024 pixel detector). The dataset consists of four 
spatial and velocity resolved images, which are amalgamated 
into one field of 89"x67"or 0.2x0.15 pc (4 1 • 10 4 x 31 ■ 10 4 
AU) given a distance to Orion of 460 pc (Ba llv et all EoOO). 
The field is centred approximately on the Becklin-Neugebauer 
(BN) object (05 h 35 m 14 s .l, -05°22'22'/9). The spatial resolution 
is 0715 (70 AU) and relative velocities between adjacent res- 
olution elements are determined to 0.3kms _1 (lcr) in bright re- 
gions to 8.4kms _1 in the weakest regions considered here. A 
detailed description of data acquisition and reduction may be 
found in IGustafsson et alJ J2003I) with further discussion in 
IGustafsson et alJ ll2006l) . We will refer to this dataset as the 

CFHT/Gri F data. 

Data in iRristensen et alJ d2006l) were also obtained at the 
CFHT, with the PUEO adaptive optics and the KIR infrared 
detector but without velocity information and covering a more 
restricted region than the CFHT/GriF data. Emission from 
the H2 1-0 S(l) line was isolated using a narrow-band fil- 
ter (/1/100). Continuum was recorded with a filter centred at 
2.183/im and subtracted from the data obtained in the S(l) 
filter , data being subje cted to standard reduction procedures 
dKristensen et alJl2006l). The dataset u sed here is that desig- 
nated as W in IRristensen et alJ d2006l) centered ~ ^"south- 
west of BN. The field is 36"x 36", the pixel scale is 07035 and 
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the resolution is 0745 (200AU). This dataset will be referred to 
as CFHT/PU EO-KIR. 

Data in Lacom be et alJ d2004l) were obtained with 
VLT/NACO. The H 2 1-0 S(l) line was isolated using a 
narrow-band fil t er, wi th continuum subtraction as described in 
lLacombe et alJ J2Q04I) . The observations used here are those 
designated as the east-south-east (ESE) field, 27"x27"with a 
pixel scale of 07027 and a resoluti on of 07 12 (55AU). T his cor- 
responds to the field analysed in IVannier et al.1 J200ll) and is 
referred to here as the VLT/NACO data. This field lies to the 
east of the CFHT/PUEO-KIR field. 

In the succeeding sections we use the new technique em- 
ploying SFs for the various datasets as follows. First we use 
brightness data extracted from the CFHT/GriF. These data 
cover the greatest physical extent of any of our datasets. At 
this stage we ignore the velocity information in these data. We 
then consider the CFHT/PUEO-KIR and VLT/NACO data. In 
the latter two cases we compare our re sults with those of the 
histogram method of Blitz & Williams) i 19971) and the Fourier 
transform method. In a third section we use the velocity data 
contained in the CFHT/GriF results. 

4. Identification of Structure in OMC1 

4.1. Use of the method with CFHT/GriF H 2 brightness 
data 

We use the CFHT/GriF observations described in Sect. [3] to 
show that the molecular gas in OMC1 contains a number of 
preferred scales, illustrating the present technique. We have 
calculated SFs and logarithmic derivatives for these data. Note 
that brightness is roughly representative of the density struc- 
ture in the excited H2. In Fig.|4]the structure function of order 5 
for the brightness and the corresponding derivatives are shown. 
Relative errors on values in the SFs are found using the law 
of propagation of errors for uncorrelated variables. Using the 
notation of Eq. Q, the variance is calculated as: 

o 2 (S p (r)) = (p 2 /N 2 ) ^(o-W)) + o-\A{r' - 5))) 

x \A(r')-A(r'-S)\ 2(p - l} (3) 

where the summation is performed over all pairs of pixels in 
the map that satisfy \6\ = r, N is the total number of such pairs 
and cr 2 (A(r')) is the variance of A(r'). Due to the large number 
of pixel pairs in our data, 2- 10 7 to 4-10 9 , the variances are very 
small. Typical values of the relative errors, cr(S P (r))/S p (r), are 
10~ 5 in the 2nd order SF, 10 4 in the 5th order SF and 10 3 
in the 10th order SF. Thus the errors are negligible in SFs of 
order up to 10, as well as in the logarithmic derivatives of these 
functions. 

Local minima in the derivatives as a function of r are 
present at r ~ 1 10AU and r ~ 1700AU, providing clear evidence 
that the brightness structure, and by implication the density dis- 
tribution, in OMC1 shows preferred scales. The positions of the 
minima indicate that the (largely) shocked regions in the field 
have two preferred sizes close to 110 and 1700AU. Two less 
pronounced minima or points of inflection are present at r ~ 
550AU and r ~ 2700AU, which suggests that clumping at these 
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Fig. 4. a) The 5th order structure function of velocity integrated 
brightness in OMC1 from GriF data from CFHT. b) Corresponding 
local logarithmic derivatives of the structure function of order 5. 



scales is also present in the field, but to a lesser degree. As 
noted earlier these sizes may be representations of diameters 
in near-circular clumps or major and minor axes in elliptically 
shaped clumps. 

In Fig. [5^ the logarithmic derivatives of SFs of order 
4,6,8,10 are displayed for the same data. This illustrates the 
tendency for the derivatives of higher order SFs to congregate 
at the same value, a property mentioned in Sect. 13 In Fig.|5J) 
the derivatives are displaced on the ordinate to give a clearer 
view of the evolution with regard to the order of the SF. As the 
order of the structure function increases, the derivatives show 
more detail. High order SFs accentuate the presence of weaker 
structure hidden in lower order SFs, resulting in more minima 
in the derivatives. Thus in the limit of very high order, so far 
as this can be achieved within the errors of observation, deriva- 
tives will essentially pick out the size of every clump in the 
region. This property is useful since the order of the structure 
function at which a structure becomes apparent is a measure of 
the prevalence of that structure. A structure that is clearly seen 
in the 2nd order SF is more dominant in the region than a struc- 
ture which is first seen in the 6th order SF, say. Examples are 
the scale sizes at r~ 550AU and 2700AU which are barely visi- 
ble in the derivatives of the 4th order SF, but are clearly seen in 
the derivatives of the 6th order SF. These are less dominant than 
the scale sizes at r~ 1 10AU, 1700AU, which are clearly appar- 
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Fig. 5. a) Local logarithmic derivatives of structure functions of order 
4,6,8,10. b) As in a) but with the graphs displaced to get a better view 
of the evolution. 

ent in the derivatives of the 4th order SF. In observational data, 
random errors in the calculated values will however increase 
when higher orders are computed, thus limiting the maximum 
value of p that may be used. Thus uncertainties in observational 
data may generally be a limiting factor in this form of data anal- 
ysis. 

4.2. Comparison with other methods using 
CFHT/PUEO-KIR and VLT/NACO data 

We first present a comparison of the method described 
here with the Fourier transform technique conduc ted on the 
CFHT/PUEO-KIR data of iKristensen etaD ll2006l) . FigureEJi 
shows the derivatives of the SFs of order 4,5,6 of that region 
indicating a preferred scale of 1600 AU. We now compare this 
with the standard two-dimensional Fourier transform technique 
commonly employed to identify scale size in images. This 
yields a power spectrum P(k) which in the case of self-similar 
structure should be a power law in k (k oc r _1 ). Deviations from 
a power law are nicely illustrated by showing k 2 P(k) vs k, as 
described in IVannier et alJ J200lb . Figure [6j> displays k 2 P(k) 
vs. scale size (k _1 ). Figure [(J) shows that there is a range of 
scales which are over-populated. This over-population peaks at 
~ 1200-1800 AU. The above results show good agreement be- 
tween the scales obtained from the power spectrum and the SF 



Fig. 6. a) Local logarithmic derivatives of structure functions of or- 
der 4,5,6 of a region South-West of BN in OMC1 calculated using 
CFHT/PUEO-KIR data, b) Fourier transform of the same set of data: 
the power spectrum multiplied by k 2 vs. the scale size. 

derivatives, at any rate in this simple case in which only a single 
scale, or a single range of scales, is present in the region. 

We now turn to a comparison with a more complicated ex- 
ample in which more than one scale is present. P referred scales 
inOMCl have previously been identified by IVannier et alJ 
J200 lh and lLacombe et all HcXM) using a variety of techniques. 
We now use the SF technique to obtain scale sizes in the 
VLT/NACO data. These may be com pared with sca l e size s 
reported in lLacombe et alJ ( 120041) and IVannier et alJ J200ll) . 
Using an area-perim e ter ana lysis and a Fourier transform tech- 
nique IVannier et alJ J200ll) found a preferred scale around 
600A U, but spanning a la rge and rather poorly determined 
range. Vannier et al. (2001) als o used the brightness histogram 
technique of Blitz & Williams d!997l) . finding a preferred scale 
of -900 AU. lLacombe et alJ J2004) with the VLT/NACO data, 
used only the histogram technique finding a preferred size of 
1100±130AU. 

The local derivatives of the 4th order SF for the VLT/NACO 
data are shown in Fig. where two scale sizes of r~500AU 
and 2000AU are evident. The lower of these scales is con- 
sistent with the findings in Vann ier et alJ J200ll) . As a further 
comparison, we created a two-dimensional Fourier transform 
to yield the power spectrum, P(k) of the same field using the 
VLT/NACO rather than the 3.6m La Silla data of lVannier et alJ 
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( 2001 ). The power spectrum wa s - as it should be - very similar 
to that reported in Vann ier et alJ J200ll) showing marked depar- 
ture from a power law in k. Data in Fig. show that the SF 
method is capable of much more precision in the identification 
of structure scale than a simple Fourier transform. 

In summary, there is reasonably good correspondence be- 
tween the standard technique of 2D Fourier transform and the 
SF method, so far as comparison can be made. Both show clear 
departure from self-similarity and a similar scale size, where 
this can be identified as in Figs. [6^ and b. There seems to be 
some discrepancy between these two technique s and the his- 
togram me thod, on the basis of w ork reported in I Vannier et alJ 
( 200 1 ) and lLacombe et alJ ( 12004 



4.3. Use of the present method with CFHT/GriF H 2 
velocity data 

We now return to the dataset covered by the CFHT/GriF obser- 
vations, as described in Sect.[3] The 5th order structure function 
of the velocity data from GriF and the corresponding deriva- 
tives are shown in Fig. [8] Because the accuracy of the veloc- 
ity in each pixel depends on the brightness in the pixel (see 
Sect. each velocity difference in Eq. Q has been weighted 
by the product of the brightness, B, in the two pixels in ques- 
tion, thus giving more weight to pixels with high brightness. 
The weighted SFs are thus 



SJr) = (B(r')B(r' - S)W) - v(r' - 5)\ p ) 



(4) 



Relative errors on the calculated values in the weighted SFs are 
typically 10~ 3 for fifth order (p=5) and therefore negligible. 

Figure [8] shows that local minima in the derivatives as a 
function of r are present at r ~ 140AU, 1500AU and 3500AU. 
This analysis therefore reveals that there is local order in the ve- 
locity field at these scales. It is interesting that these scales are 
respectively disk scales (140 AU) and outflow scales (>1000 
AU), the latter being associated with (re-)injection of energy 
into the system, that is, a driving scale. These scales are ap- 
proximately the same as those identified in the brightness data, 
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Fig. 8. a) The 5th order structure function of velocities in OMC1 from 
CFHT/GriF. b) Corresponding local logarithmic derivatives of the 5th 
order structure function. 



which were found to be at 110 and 1700 AU, although there 
appears to be an additional larger scale at 3500 AU in velocity. 
The weak structure in the brightness at i~ 550AU and 2700AU 
is however not seen in the velocity. An analysis of the veloc- 
ity data has also been performed without brightness weighting 
with essentially the same result as in Fig.|H] 

4.4. Systematic errors associated with estimates of 
structure scale 

There are systematic errors inherent in the present method. As 
already noted, derivatives of higher order SFs tend to pick out 
scales which are too small by as much as 10-20%. It turns out 
that when we consider a distribution of clump sizes about some 
mean value, this causes an overestimation of the mean clump 
size as we discuss below. These two effects tend to cancel but 
to an extent which cannot be quantified. 

The detection of a preferred scale size implies that there 
is a range of sizes around the derived value. In common with 
other methods, our method does not give any measure of the 
extent of such a range. Here we test the effect of introducing a 
distribution of clump sizes about some mean preferred scale us- 
ing artificial images. In order to illustrate the behaviour of such 
systems, we assume that the sizes of clumps are gaussian dis- 
tributed about the mean. Using simulations similar to those in 
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Sect. [2]we have created a number of maps with 200 randomly 
positioned clumps. Three mean sizes of 40, 60 and 80 pixels 
were used. We have found that when the half-width is smaller 
than ~5% of the mean value, the SFs are indistinguishable from 
the original of only one clump size and the minimum of the 
SF derivatives is found at the mean value. Increasing the half- 
width beyond ~5% causes the minimum to move towards larger 
scales. If the half-width of the distribution is 10% (20%) 25% 
of the mean value the minimum is found at 1.07 (1.25) 1.50 
times the mean value. For example, if the data used in Fig. [6] 
consist of clumps distributed as a gaussian with a halfwidth 
of 20% of the mean value, the detected value of 1600AU cor- 
responds to a true mean value of the sizes of structures of ~ 
1300AU. This analysis only conveys the sensitivity of the sys- 
tematic error to the width of the distribution but does not allow 
an estimation of the width from observations. The distribution 
of clump sizes about some preferred scale is directly linked to 
the IMF and it would be useful to develop methods which were 
capable of making a reliable width estimate. 

5. Conclusion 

In conclusion the technique presented here, whilst suffering 
from some degree of uncertainty and imprecision to an extent 
common to all methods, has proved valuable for analysis of 
images in our own work. With the present and future advent of 
higher spatial resolution maps in the radio region of the spec- 
trum, probing the star-forming scale, for example obtained with 
the Plateau de Bure interferometer, the Submillimeter Array, 
CARMA and ALMA, this tool may prove useful to an increas- 
ing community of observers. 
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